Probing the Low-Energy Electronic Structure of Complex Systems by ARPES 
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Angle-resolved photoemission spectroscopy (ARPES) is one of the most direct methods of study- 
ing the electronic structure of solids. By measuring the kinetic energy and angular distribution of the 
electrons photoemitted from a sample illuminated with sufficiently high-energy radiation, one can 
gain information on both the energy and momentum of the electrons propagating inside a material. 
This is of vital importance in elucidating the connection between electronic, magnetic, and chemical 
structure of solids, in particular for those complex systems which cannot be appropriately described 
within the independent-particle picture. The last decade witnessed significant progress in this tech- 
nique and its applications, thus ushering in a new era in photoelectron spectroscopy; today, ARPES 
experiments with 2 meV energy resolution and 0.2° angular resolution are a reality even for photoe- 
mission on solids. In this paper we will review the fundamentals of the technique and present some 
illustrative experimental results; we will show how ARPES can probe the momentum-dependent 
electronic structure of solids providing detailed information on band dispersion and Fermi surface, 
as well as on the strength and nature of those many-body correlations which may profoundly affect 
the one-electron excitation spectrum and, in turn, determine the macroscopic physical properties. 



I. INTRODUCTION 

Photoelectron spectroscopy is a general term that 
refers to all those techniques based on the application 
of the photoelectric effect originally observed by Hertz 
and later explained as a manifestation of the quan- 
tum nature of light by Einstein 2] , who recognized that 
when light is incident on a sample an electron can ab- 
sorb a photon and escape from the material with a max- 
imum kinetic energy Ekin = hv—cf) (where v is the photon 
frequency and 0, the material work function, is a mea- 
sure of the potential barrier at the surface that prevents 
the valence electrons from escaping, and is typically 4- 
5 eV in metals). In the following, we will show how 
the photoelectric effect also provides us with deep in- 
sights into the quantum description of the solid state. 
In particular, we will give a general overview of angle- 
resolved photoemission spectroscopy (so-called ARPES), 
a highly advanced spectroscopic method that allows the 
direct experimental study of the momentum-dependent 
electronic band structure of solids. For a further discus- 
sion of ARPES and other spectroscopic techniques based 
on the detection of photoemitted electrons, we refer the 
reade^^^^^tensiw^^^^^^jglaU^OT^t^j^j^ 

As we will see in detail throughout the paper and in 
particular in Sec. lIIII due to the complexity of the pho- 
toemission process in solids the quantitative analysis of 
the experimental data is often performed under the as- 
sumption of the independent-particle picture and of the 
sudden approximation (i.e., disregarding the many-body 
interactions as well as the relaxation of the system dur- 
ing the photoemission itself). The problem is further 
simplified within the so-called three-step model (Fig.^), 
in which the photoemission event is decomposed in three 
independent steps: optical excitation between the initial 
and final &itZfcBloch eigenstates, travel oi the excited elec- 



tron to the surface, and escape of the photoelectron into 
vacuum after transmission through the surface potential 
barrier. This is the most common approach, in partic- 
ular when photoemission spectroscopy is used as a tool 
to map the electronic band structure of solids. However, 
from the quantum-mechanical point of view photoemis- 
sion should not be described in terms of several indepen- 
dent events but rather as a one-step process (Fig.^): 
in terms of an optical transition (with probability given 
by Ea. ll2() between initial and final states consisting of 
many-body wave functions that obey appropriate bound- 
ary conditions at the surface of the solid. In particular 
(see Fig.|21l, the initial state should be one of the possible 
A^-electron eigenstates of the semi-infinite crystal, and 
the final state must be one of the eigenstates of the ion- 
ized (TV— l)-electron semi-infinite crystal; the latter has 
also to include a component consisting of a propagat- 
ing plane-wave in vacuum (to account for the escaping 
photoelectron) with a finite amplitude inside the crystal 
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FIG. 1: Pictorial representation of three-step and one-step 
model description of the photoemission process (from Ref. ilTi) . 
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(to provide some overlap with the initial state). Fur- 
thermore, as expressed by Eg . 1121 which does represent a 
complete one-step description of the problem, in order 
for an electron to be photoemitted in vacuum not only 
there must be a finite overlap between the amplitude of 
initial and final states, but also the following energy and 
momentum conservation laws for the impinging photon 
and the A^-electron system as a whole must be obeyed: 

Ef-E^^hu (1) 
kf-kf = k„, (2) 

where the indexes i and / refer to initial and final state, 
respectively, and k/j^ is the momentum of the incoming 
photon. Note that, in the following, in proceeding with 
the more detailed analysis of the photoemission process 
as well as its application to the study of the momentum- 
dependent electronic structure of solids (in terms of both 
conventional band mapping as well as many-body ef- 
fects), we will mainly restrict ourselves to the context 
of the three-step model and the sudden approximation. 

II. KINEMATICS OF PHOTOEMISSION 

The energetics and kinematics of the photoemission 
process are shown in Fig|31 andQ] while the geometry of 
an ARPES experiment is sketched in Fig|S^. A beam 
of monochromatized radiation supplied either by a gas- 
discharge lamp or by a synchrotron beamline is incident 
on a sample (which has to be a properly aligned sin- 
gle crystal in order to perform angle or, equivalently, 
momentum-resolved measurements). As a result, elec- 
trons are emitted by photoelectric effect and escape in 
vacuum in all directions. By collecting the photoelec- 
trons with an electron energy analyzer characterized by 
a finite acceptance angle, one measures their kinetic en- 
ergy Ekin for a given emission direction. This way, the 
wave vector or momentum K = p/fi of the photoelectrons 
in vacuum is also completely determined: its modulus is 
given hy K — y/2/rnEki^/h and its components parallel 
(K|| = Kx + Ky) and perpendicular (K_l = K^) to the 
sample surface are obtained in terms of the polar (jS) and 
azimuthal (ip) emission angles defined by the experiment: 



= - ^ 2mEkin sin i? cos (p (3) 

h 

Ky = i a/ 2mEkin sin -& sin ip (4) 

= \ \J 2mEkin cos (5) 

The goal is then to deduce the electronic dispersion re- 
lations -E(k) for the solid left behind, i.e. the relation 
between binding energy Eb and momentum k for the 
electrons propagating inside the solid, starting from Ekin 
and K measured for the photoelectrons in vacuum. In 
order to do that, one has to exploit the total energy and 
momentum conservation laws (Eq.^andl^l respectively). 
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FIG. 2: Initial (left) and final (right) eigenstates for the semi- 
infinite crystal. Left: (a) surface resonance; (b) surface Shock- 
ley state situated in a gap of the bulk band structure; (c) bulk 
Bloch state. Right; (d) surface resonance; (e) in gap evanes- 
cent state; (f) bulk Bloch final state (from Ref. l35ii) . 

Within the non-interacting electron picture, it is par- 
ticularly straightforward to take advantage of the energy 
conservation law and relate, as pictorially described in 
Fig.El the kinetic energy of the photoelectron to the bind- 
ing energy Eb of the electronic-state inside the solid: 

Ek^n ^ hl^ ~ - \Eb\ (6) 

More complex, as we will discuss below, is to gain full 
knowledge of the crystal electronic momentum k. Note, 
however, that the photon momentum can be neglected 
in Eq.|2] at the low photon energies most often used in 
ARPES experiments {hi^K 100 eV), as it is much smaller 
than the typical Brillouin-zone dimension 27r/a of a solid 
(see Sec. lVIl for more details). Thus, as shown in Fig.^ 
within the three-step model description (see also Sec. llllll . 
the optical transition between the bulk initial and final 
states can be described by a vertical transition in the 
reduced-zone scheme (kj— = 0), or equivalently by a 
transition between momentum-space points connected by 
a reciprocal-lattice vector G in the extended-zone scheme 
(k/ — ki = G). In regard to Eq.Q] ancl[21 and the deeper 
meaning of the reciprocal-lattice vector G note that, as 
emphasized by Mahan in his seminal paper on the the- 
ory of photoemission in simple metals 36], "m a nearly- 
free- electron gas, optical absorption may be viewed as a 
two-step process. The absorption of the photon provides 
the electron with the additional energy it needs to get to 
the excited state. The crystal potential imparts to the 
electron the additional momentum it needs to reach the 
excited state. This momentum comes in multiples of the 
reciprocal-lattice vectors G. So in a reduced zone picture, 
the transitions are vertical in wave-vector space. But in 
photoemission, it is more useful to think in an extended- 
zone scheme." On the contrary in an infinite crystal with 
no periodic potential (i.e., a truly free-electron gas sce- 
nario lacking of any periodic momentum structure), no 
k-conserving transition is possible in the limit k/i^ = 0, 
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FIG. 3: Energetics of the photoemission process (from 
Ref. H^ . The electron energy distribution produced by the in- 
coming photons, and measured as a function of the kinetic en- 
ergy Ekin of the photoelectrons (right), is more conveniently 
expressed in terms of the binding energy Eb (left) when one 
refers to the density of states in the solid {Eo,b = at Ep). 



as one cannot go from an initial to a final state along 
the same unperturbed free-electron parabola without an 
external source of momentum. In other words, direct 
transitions are prevented because of the lack of appropri- 
ate final states (as opposed to the periodic case of Fig.^. 
Then again the problem would be quite different if the 
surface was more realistically taken into account, as in a 
one-step model description of a semi-infinite crystal. In 
fact, while the surface does not perturb the translational 
symmetry in the x-y plane and k|| is conserved to within 
a reciprocal lattice vector G||, due to the abrupt poten- 
tial change along the z axis the perpendicular momentum 
kj_ is not conserved across the sample surface (i.e., k_L is 
not a good quantum number except than deeply into the 
solid, contrary to k| | ) . Thus, the surface can play a direct 
role in momentum conservation, delivering the necessary 
momentum for indirect transitions even in absence of the 
crystal potential (i.e., the so-called surface photoelectric 
effect] see also Ea.lT^and the related discussion). 

Reverting to the three-step model direct-transition de- 
scription of Fig.21 the transmission through the sample 
surface is obtained by matching the bulk Bloch eigen- 
states inside the sample to free-electron plane waves in 
vacuum. Because of the translational symmetry in the 
x-y plane across the surface, from these matching condi- 
tions it follows that the parallel component of the elec- 
tron momentum is actually conserved in the process: 
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k|| = K|| = -^j2mEkin ■ sini? 



(7) 



electron crystal momentum in the extended-zone scheme 
(upon going to larger angles, one actually probes elec- 
trons with k|| lying in higher-order Brillouin zones; by 
subtracting the corresponding reciprocal-lattice vector 
G||, the reduced electron crystal momentum in the first 
Brillouin zone is obtained). As for the determination 
of k_L, which is not conserved but is also needed in or- 
der to map the electronic dispersion ^^(k) vs the total 
crystal wave vector k, a different approach is required. 
As a matter of fact, several specific experimental meth- 
ods for absolute three dimensional band mapping have 
been developed 0, [sl, [s^, which however are rather 
complex and require additional and/or complementary 
experimental data. Alternatively, the value of kj^ can be 
determined if some a priori assumption is made for the 
dispersion of the electron final states involved in the pho- 
toemission process; in particular, one can either use the 
results of band structure calculations, or adopt a nearly- 
free-electron description for the final bulk Bloch states: 



2m 



(8) 



where once again the electron momenta are defined in 
the extended-zone scheme, and Eq corresponds to the 
bottom of the valence band as indicated in Fig.^ (note 
that both Eq and Ef are referenced to the Fermi energy 
Ep, while is referenced to the vacuum level Ey). 

Because E f — E}^in + 4' ^nd ?i^k|/2m = iJfei„ sin^ i9, which 
follow from Fig. 0] and Eq.[7| one obtains from Eq.|SJ 
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k± = - V 2m{Eki„ cos2 ■& + ¥()) 



(9) 



Here Vq = \Eo\-\-4> is the inner potential, which corre- 
sponds to the energy of the bottom of the valence band 
referenced to vacuum level Ey . From Eq.EI and the mea- 
sured values of E^i^ and i?, if Vq is also known, one can 
then obtain the corresponding value of kj_ . As for the de- 
termination of Vq, three methods are generally used: (i) 
optimize the agreement between theoretical and experi- 
mental band mapping for the occupied electronic state; 




where k|| is the component parallel to the surface of the 
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FIG. 4: Kinematics of the photoemission process within the 
three-step nearly-free-electron final state model: (a) direct op- 
tical transition in the solid (the lattice supplies the required 
momentum); (b) free-electron final state in vacuum; (c) cor- 
responding photoelectron spectrum, with a background due 
to the scattered electrons (-Bo,i,/=0 at Ep)- From Ref.ls^ 
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FIG. 5; Normal and grazing emission ARPES spectra from 
Ag(lOO) measured with photon energies specifically chosen to 
give rise to peaks with the same binding energy (from Ref . liol) . 



(ii) set Vo equal to the theoretical zero of the muffin tin 
potential used in band structure calculations; (iii) infer 
Vo from the experimentally observed periodicity of the 
dispersion E{'k±). The latter is actually the most conve- 
nient method as the experiment can be realized by simply 
detecting the photoelectrons emitted along the surface 
normal (i.e., K|| =0) while varying the incident photon 
energy and, in turn, the energy Ekin of the photoelec- 
trons and thus (see Eq.lSJ. Note that the nearly-free 
electron approximation for the final states is expected to 
work well for materials in which the Fermi surface has 
a simple spherical (free-electron-like) topology such as in 
the alkali metals, and for high-energy final states in which 
case the crystal potential is a small perturbation (even- 
tually the final-state bands become so closely spaced in 
energy to form a continuum, and the details of the fi- 
nal states become unimportant). However this approxi- 
mation is often used also for more complicated systems, 
even if the initial states are not free electron-like. 

A particular case in which the uncertainty in is less 
relevant is that of the low-dimensional systems character- 
ized by an anisotropic electronic structure and, in par- 
ticular, a negligible dispersion along the z axis (i.e., the 
surface normal, see Fig.|^). The electronic dispersion is 
then almost exclusively determined by k|| (as in the case 
of many transition metal oxides, such as for example the 
two-dimensional copper oxide superconductors As 
a result, one can map out in detail the electronic dis- 
persion relations E(]s.) simply by tracking, as a function 
of K||, the energy position of the peaks detected in the 
ARPES spectra for different take-off angles (as in Fig.EjD, 
where both direct and inverse photoemission spectra for 
a single band dispersing through the Fermi energy Ep 



are shown). Furthermore, as an additional bonus of the 
lack of z dispersion, one can directly identify the width 
of the photoemission peaks with the lifetime of the pho- 
tohole '4l"|, which contains information on the intrinsic 
correlation effects of the system and is formally described 
by the imaginary part of the electron self energy (see 
Sec. lIV|l . On the contrary, in 3D systems the linewidth 
contains contributions from both photohole and photo- 
electron lifetimes, with the latter reflecting final state 
scattering processes and thus the finite probing depth; 
as a consequence, isolating the intrinsic many-body ef- 
fects becomes a much more complicated problem. 

What just discussed for the lifetime can be easily seen 
from the expression for the FWHM of an ARPES line- 
shape for a single nearly- free electron-like band |4l| : 

r.- , r. 
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(10) 



Here F f and F; are the inverse lifetime of photoelectron 
and photohole in the final and initial states, respectively, 
and Vi and Vf are the corresponding group velocities (e.g., 
hvi± = dEi/dk±). Note in particular that: (i) for initial 
states very close to Ep, F^ — >0 and the linewidth reflects 
only the lifetime of the final state F/; (ii) Ea.llOlsimplifies 
considerably in the case of a material characterized by a 
two dimensional electronic structure, for which |i;i^|~0; 
as a result, the final-state lifetime contribution vanishes: 
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(11) 



Furthermore, depending on the sign of Vi\\, the measured 
linewidth can be compressed or expanded with respect to 
the intrinsic value of the inverse lifetime F^ . The two lim- 
iting cases mentioned above are beautifully exemplified 
by the data from the three-dimensional system Ag(lOO) 
presented in Fig.|Sl ^3 ■ While the normal incidence spec- 
trum is dominated by Fy ^ Fj and is extremely broad, 
the grazing incidence data from a momentum space re- 
gion characterized by Vi± = 0, < and large, and 
/c|| small (which result in a compression factor C = 0.5), 
exhibit a linewidth which is even narrower than the in- 
trinsic inverse lifetime F^. Note that this does not im- 
ply any fundamental violation of the basic principles of 
quantum mechanics, but is just a direct consequence of 
the kinematics constrains of the photoemission process. 



III. THREE-STEP MODEL AND SUDDEN 
APPROXIMATION 

To develop a formal description of the photoemission 
process, one has to calculate the transition probabil- 
ity Wfi for an optical excitation between the A'^-electron 
ground state \E'f and one of the possible final states 
This can be approximated by Fermi's golden rule: 

-""'l^SiEf -E^ -hi.) (12) 



Wfi 
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where =E^-^-E^ and Ef = E^-^+Ekin are the 
initial and final-state energies of the A^-particle system 
{E^ is the binding energy of the photoelectron with ki- 
netic energy Ekin and momentum k). The interaction 
with the photon is treated as a perturbation given by: 

= 7r^(A-p + p-A) = — A-p (13) 
Imc mc 

where p is the electronic momentum operator and A 
is the electromagnetic vector potential (note that the 
gauge $ = was chosen for the scalar potential $, 
and the quadratic term in A was dropped because in 
the linear optical regime it is typically negligible with 
respect to the linear terms) . In Eg. 1131 we also made 
use of the commutator relation [p, A] = —ihV ■ A and 
dipole approximation [i.e., A constant over atomic di- 
mensions and therefore V • A — 0, which holds in the 
ultraviolet]. Although this is a routinely used approxi- 
mation, it should be noted that V-A might become im- 
portant at the surface where the electromagnetic fields 
may have a strong spatial dependence. This surface 
photoemission contribution, which is proportional to 
(e — 1) where e is the medium dielectric function, can 
interfere with the bulk contribution resulting in asym- 
metric lineshapes for the bulk direct-transition peaks 
E3) 13 • At this point, a more rigorous approach 
is to proceed with the so-called one-step model (Fig.^), 
in which photon absorption, electron removal, and elec- 
tron detection are treated as a single coherent process 

[illillillillilliiiSElilillllElS 

In this case bulk, surface, and vacuum have to be in- 
cluded in the Hamiltonian describing the crystal, which 
implies that not only bulk states have to be considered 
but also surface and evanescent states, and surface res- 
onances (see Fig. 121. Note that, under the assumption 
V-A = 0, from Eg. 1131 and the commutation relation 
[7Jo,p] = ifi\/V (where Hq = p'^/2ni + V is the unper- 
turbed Hamiltonian of the semi-infinite crystal) it follows 
that the matrix elements appearing in Eg . 1121 are propor- 
tional to (f^lA-VFl^f ). This shows explicitly that for 
a true free-electron like system it would be impossible 
to satisfy simultaneously energy and momentum conser- 
vation laws inside the material because there W = 0. 
The only region where electrons could be photoexcited 
is at the surface where dV/dz ^ 0, which gives rise to 
the so-called surface photoelectric effect. However, due 
to the complexity of the one-step model, photoemission 
data are usually discussed within the three-step model 
(Fig.^), which although purely phenomenological has 
proven to be rather successful 49^ ^SO^ ,6^ . Within this 
approach, the photoemission process is subdivided into 
three independent and sequential steps: 

(i) Optical excitation of the electron in the hulk. 

(ii) Travel of the excited electron to the surface. 

(iii) Escape of the photoelectron into vacuum. 



The total photoemission intensity is then given by the 
product of three independent terms: the total probabil- 
ity for the optical transition, the scattering probability 
for the travelling electrons, and the transmission proba- 
bility through the surface potential barrier. Step (i) con- 
tains all the information about the intrinsic electronic 
structure of the material and will be discussed in detail 
below. Step (ii) can be described in terms of an effective 
mean free path, proportional to the probability that the 
excited electron will reach the surface without scattering 
(i.e, with no change in energy and momentum). The in- 
elastic scattering processes, which determine the surface 
sensitivity of photoemission (see Sec. lVI|) . give rise to a 
continuous background in the spectra which is usually ig- 
nored or subtracted. Step (iii) is described by a transmis- 
sion probability through the surface, which depends on 
the energy of the excited electron and the material work 
function (j) (in order to have any finite escape probability 
the condition fi,^k^/2TO> |£'o|+(/) must be satisfied). 

In evaluating step (i), and therefore the photoemis- 
sion intensity in terms of the transition probability Wfi, 
it would be convenient to factorize the wavefunctions in 
Eq.Elinto photoelectron and (A^— l)-electron terms, as 
we have done for the corresponding energies. This how- 
ever is far from trivial because during the photoemission 
process itself the system will relax. The problem sim- 
plifies within the sudden approximation^ which is exten- 
sively used in many-body calculations of the photoemis- 
sion spectra from interacting electron systems, and is in 
principle applicable only to high kinetic- energy electrons. 
In this limit, the photoemission process is assumed to be 
sudden, with no post-coUisional interaction between the 
photoelectron and the system left behind (in other words, 
an electron is instantaneously removed and the effective 
potential of the system changes discontinuously at that 
instant). The final state can then be written as: 

■^'J =Ac^)-^^-^ (14) 

where A is an antisymmetric operator that properly an- 
tisymmetrizes the iV-electron wavefunction so that the 
Pauli principle is satisfied, is the wavefunction of the 
photoelectron with momentum k, and is the fi- 

nal state wavefunction of the (A^— l)-electron system left 
behind, which can be chosen as an excited state with 
eigenfunction 5"^"^ and energy E^~^. The total transi- 
tion probability is then given by the sum over all possible 
excited states m. Note, however, that the sudden approx- 
imation is inappropriate for low kinetic energy photoelec- 
trons, which may need longer than the system response 
time to escape into vacuum. In this case, the so-called 
adiabatic limit, one can no longer factorize in two 
independent parts and the detailed screening of photo- 
electron and photohole has to be taken into account [g^ ■ 
For the initial state, let us first assume for simplicity 
that ^'f is a single Slater determinant (i.e., Hartree-Fock 
formalism), so that we can write it as the product of a 
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FIG. 6: (a) Geometry of an ARPES experiment; the emission direction of the photoelectron is specified by the polar (i?) 
and azimuthal (ip) angles. Momentum resolved one-electron removal and addition spectra for: (b) a non-interacting electron 
system (with a single energy band dispersing across the Fermi level); (c) an interacting Fermi liquid system. The corresponding 
ground-state (T = K) momentum distribution function n(k) is also shown, (c) Bottom right: photoelectron spectrum of 
gaseous hydrogen and the ARPES spectrum of solid hydrogen developed from the gaseous one (from Ref.|^ andJSa). 



one-electron orbital (f)f and an (A^— l)-particle term: 



N 



(15) 



More generally, however, ^ should be expressed as 
^ff-i=Ck4'f , where Ck is the annihilation operator for 
an electron with momentum k. This also shows that 
is not an eigenstate of the (A^— 1) particle Hamil- 
tonian, but is just what remains of the A^-particle wave- 
function after having pulled out one electron. At this 
point, we can write the matrix elements in Eg . 1121 as: 



(16) 



where ((/)^|i?i„f |0^) = M^j is the one-electron dipole ma- 
trix element, and the second term is the (A^— l)-electron 
overlap integral. Here, we replaced ^'^ ~^ with an eigen- 
state SIS discussed above. The total photoemission 
intensity measured as a function of Ekin at a momentum 
k, namely /(k, Ekin) — X)/ i''^f,i^ then proportional to: 



where |c„ 



\''5{Ek^n+E, 



N-1 



E^^hiy) (17) 



7V-l|^W-lN|2 



is the probability that 



the removal of an electron from state i will leave the 
(A^— l)-particle system in the excited state to. From here 



we see that, if '^^'^ = ^mo^ '^^'^ particular TO = mo, 
the corresponding Icmg^jp will be unity and all the others 
Cm,i zero; in this case, if also A/^ ^t^O, the ARPES spec- 
tra will be given by a delta function at the Hartree-Fock 
orbital energy E^ = — Ck, as shown in Fig.|S|D (i.e., non- 
interacting particle picture). In the strongly correlated 
systems, however, many of the |cm,ip will be different 
from zero because the removal of the photoelectron re- 
sults in a strong change of the system effective potential 
and, in turn, will have an overlap with many of 

the eigenstates Therefore, the ARPES spectra 

will not consist of single delta functions but will show a 



main line and several satellites according to the number 
of excited states to created in the process (Fig.|B|:) . 

What discussed above is very similar to the situation 
encountered in photoemission from molecular hydrogen 
|63| in which not simply a single peak but many lines 
separated by few tenths of eV from each other are ob- 
served (solid line in Fig.E):, bottom right). These so- 
called 'shake-up' peaks correspond to the excitations of 
the different vibrational states of the molecule. In 
the case of solid hydrogen (dashed hne in Fig.lHt, bot- 
tom right), as discussed by [6ll |. the vibrational excita- 
tions would develop in a broad continuum while a sharp 
peak would be observed for the fundamental transition 
(from the ground state of the II2 to the one of the 
molecule). Note that the fundamental line would also be 
the only one detected in the adiabatic limit, in which case 
the (A^— l)-particle system is left in its ground state. 



IV. ONE-PARTICLE SPECTRAL FUNCTION 

In the discussion of photoemission on solids, and in 
particular on the correlated electron systems in which 
many |Cm,iP in Eg. 1171 are different from zero, the most 
powerful and commonly used approach is based on the 
Green's function formahsm jH El S Elii El . In 
this context, the propagation of a single electron in a 
many-body system is described by the time-ordered one- 
electron Green's function Q{t — t'), which can be in- 
terpreted as the probability amplitude that an electron 
added to the system in a Bloch state with momentum 
k at a time zero will still be in the same state after a 
time |t— By taking the Fourier transform, G{t—t') can 
be expressed in energy-momentum representation result- 
ing in g{k,uj) = G+{k,uj) + G-{k,uj), where G+{k,uj) 
and G~ (k, lo) are the one-electron addition and removal 
Green's function, respectively. At r = 0: 



G±(k,cj) 



N±l\ 



I*, 



N\\2 



E, 



N±l 



Ef ± iri 



(18) 
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where the operator cj^ 



^ko-) creates (anni- 



hilates) an electron with energy uj, momentum k, and 
spin cr in the A''-particle initial state ^f, the summation 
runs over all possible (A^±l)-particle eigenstates ^^^^ 
with eigenvalues E^^^ , and 77 is a positive infinitesi- 
mal (note also that fi'om here on wc will take h = 1). 
In the limit 77 — > 0+ one can make use of the identity 
(xzLir])~^ = V{l/x)^iTr6{x), where V denotes the prin- 
ciple value, to obtain the one-particle spectral function 
A{k, uj) = A+ (k, uj)+A~ (k, oj) = -(l/7r)Im G(k, w), with: 

A±(k,c.)=^|(*^±i|c±|vI/f)PJ(..-i?^±Vi?f) (19) 

m 

and G{k,Lu) = G+(k, w) -t- [G" (k, w)]*, which defines 
the retarded Green's function. Note that A~(k, w) and 
A+(k, w) define the one-electron removal and addition 
spectra which one can probe with direct and inverse pho- 
toemission, respectively. This is evidenced, for the di- 
rect case, by the comparison between the expression for 
A~ (k, oj) and Eq.^lfor the photoemission intensity (note 
that in the latter ^f^"^ ^Ck^'f^ and the energetics of the 
photoemission process has been explicitly accounted for). 
Finite temperatures effect can be taken into account by 
extending the Green's function formalism just introduced 
to T^O (see, e.g., Ref. ls^ . In the latter case, by invok- 
ing once again the sudden approximation the intensity 
measured in an ARPES experiment on a 2D single-band 
system can be conveniently written as: 



/(k,c^) =/o(k, j,,A)f{uj)A{-k,iu) 



(20) 



where k = ky is the in-plane electron momentum, uj is 
the electron energy with respect to the Fermi level, and 
Jo(k, 7^, A) is proportional to the squared one-electron 
matrix element |My and therefore depends on the elec- 
tron momentum, and on the energy and polarization of 
the incoming photon. We also introduced the Fermi func- 
tion f{uj) = (e"/'^^"^-!-!)""'^ which accounts for the fact that 
direct photoemission probes only the occupied electronic 
states. Note that in Eq.EOlwe neglected the presence of 
any extrinsic background and the broadening due to the 
energy and momentum resolution, which however have to 
be carefully considered when performing a quantitative 
analysis of the ARPES spectra (see Sec.lVland Eg .1771). 

The corrections to the Green's function originating 
from electron-electron correlations can be conveniently 
expressed in terms of the electron proper self energy 
E(k, Lu) — S'(k, CL')-HS"(k, cu). Its real and imaginary part 
contain all the information on the energy renormaliza- 
tion and lifetime, respectively, of an electron with band 
energy e/c and momentum k propagating in a many-body 
system. The Green's and spectral functions expressed in 
terms of the self energy are then given by: 



1 



Aik,io) 



- fik - S(k,u;) 
S"(k,c.) 



[c^-ek-S'(k,c^)]2 + p"(k,c^)]2 



(21) 
(22) 



Because G{t,t') is a linear response function to an ex- 
ternal perturbation, the real and imaginary parts of its 
Fourier transform G(k, lu) have to satisfy causality and, 
therefore, are related by Kramers-Kronig relations. This 
implies that if the full A(k,uj) = — (l/7r)IniG(k, w) is 
available from photoemission and inverse photoemission, 
one can calculate Re G(k, oj) and then obtain both the 
real and imaginary parts of the self energy directly from 
Eq.E] However, due to the lack of high-quality inverse 
photoemission data, this analysis is usually performed 
only using ARPES spectra by taking advantage of cer- 
tain approximations (such as, e.g., particle-hole symme- 
try within a narrow energy range about Ep \7(^). 

In general, the exact calculation of E(k, w) and, in 
turn, of A(k, u) is an extremely difficult task. In the 
following, as an example we will briefly consider the in- 
teracting FL case [t^, IZ^- Let us start from the 
trivial I](k, w) =0 non-interacting case. The A^-particle 
eigenfunction is a single Slater determinant and we 
always end up in a single eigenstate when removing 
or adding an electron with momentum k. Therefore, 
G(k, uj) — l/{uj—ekHr]) has only one pole for each k, and 
A(k, uj) = S{uj—ek) consists of a single line at the band en- 
ergy efe (as shown in Fig.O)). In this case, the occupation 
numbers rij^^ — cl.^Cy^^ are good quantum numbers and 
for a metallic system the momentum distribution [i.e., 
the expectation value ri,(k) = (riko-), quite generally in- 
dependent of the spin a for nonmagnetic systems], is 
characterized by a sudden drop from 1 to at h — kp 
(Fig. 03, top), which defines a sharp Fermi surface. If we 
now switch on the electron-electron correlation adiabat- 
ically, (so that the system remains at equilibrium), any 
particle added into a Bloch state has a certain probabil- 
ity of being scattered out of it by a collision with an- 
other electron, leaving the system in an excited state in 
which additional electron-hole pairs have been created. 
The momentum distribution n(k) will now show a dis- 
continuity smaller than 1 at kj? and a finite occupation 
probability for k > kj? even at T = (Fig.lHJ:, top). As 
long as n(k) shows a finite discontinuity Zk > at k = kj? , 
we can describe the correlated Fermi sea in terms of well 
defined quasiparticles, i.e. electrons dressed with a man- 
ifold of excited states, which are characterized by a pole 
structure similar to the one of the non-interacting sys- 
tem but with renormalized energy Ck and mass m*, and 
a finite lifetime Tk = l/Fk- In other words, the proper- 
ties of a FL are similar to those of a free electron gas 
with damped quasiparticles. As the bare-electron char- 
acter of the quasiparticle or pole strength (also called 
coherence factor) is .^k < 1 and the total spectral weight 
must be conserved (see Eg . 125(1 . we can separate G(k, w) 
and A(k, uj) into a coherent pole part and an incoherent 
smooth part without poles 

G(k,w 
A{k,uj)^Zk 



: po i 
III: 



^k 



UJ - ek + zFk 
(w-ek)2+ri 



Ginch 



(23) 
(24) 
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FIG. 7: (a) Mirror plane emission from a d^2_y2 orbital, (b) 
Sketch of the optical transition between atomic orbitals with 
different angular momenta (the harmonic oscillator wavefunc- 
tions are here used for simplicity) and free electron wave- 
functions with different kinetic energies (from Ref.[l3). (c) 
Calculated photon energy dependence of the photoionization 
cross-sections for Cu 3d and O 2p atomic levels (from Ref. l?^ . 



where Zk = (l-ff)-i, 



and the self energy and its derivatives are evaluated at 
uj = eit_. It should be emphasized that the FL descrip- 
tion is valid only in proximity to the Fermi surface and 
rests on the condition ek — for small (w — /x) and 
(k-ki^). Furthermore, Fk (x [(TrfcRp^+Jgw-Ml^l for a FL 
system in two or more dimensions [iJ, [Z3 , although ad- 
ditional logarithmic corrections should be included in the 
two-dimensional case [t^ . By comparing the electron re- 
moval and addition spectra for a FL of quasiparticles with 
those of a non-interacting electron system (in the lattice 
periodic potential), the effect of the self-energy correc- 
tion becomes evident (see Fig.O; and E|d, respectively). 
The quasiparticle peak has now a finite lifetime (due to 
E"), and it sharpens up rapidly thus emerging from the 
broad incoherent component upon approaching the Fermi 
level, where the lifetime is infinite corresponding to a well 
defined quasiparticle [note that the coherent and incoher- 
ent part of A(k, Lu) represent the main line and satellite 
structure discussed in the previous section and shown in 
Fig.Et, bottom right]. Furthermore, the peak position 
is shifted with respect to the bare band energy ek (due 
to S'): as the quasiparticle mass is larger than the band 
mass because of the dressing (m* > m) , the total disper- 
sion (or bandwidth) will be smaller (|£k| < |ek|)- 

Among the general properties of the spectral function 
there are also several sum rules. A fundamental one, 
which in discussing the FL model was implicitly used to 
state that / dojAch = -^k and J dujAinch = 1 — ^k (where 
Ach and Ainch refer to coherent and incoherent parts of 



the spectral function, respectively), is the following: 



dujA{k,uj) = l (25) 



which reminds us that A(k, uj) describes the probability 
of removing/adding an electron with momentum k and 
energy w to a many-body system. However, as it also 
requires the knowledge of the electron addition part of 
the spectral function, it is not so useful in the analysis of 
ARPES data. A sum rule more relevant to this task is: 



du!f{uj)A{k,uj) = n{k) 



(26) 



which solely relates the one-electron removal spectrum to 
the momentum distribution n(k). When electronic cor- 
relations are important and the occupation numbers are 
no longer good quantum numbers, the discontinuity at 
ki? is reduced (as discussed for the FL case) but a drop 
in n(k) is usually still observable even for strong corre- 
lations [t^. By tracking in fc-space the loci of steepest 
descent of the experimentally determined n(k), i.e. max- 
ima in |Vk 'T-(k)|, one may thus identify the Fermi surface 
even in those correlated systems exhibiting particularly 
complex ARPES features. However, great care is neces- 
sary in making use of Eq.l^Slbecause the integral of Elq.l2UI 
does not give just ?i(k) but rather /o(k, i/, A)n(k) 



V. MATRIX ELEMENTS AND FINITE 
RESOLUTION EFFECTS 

As discussed in the previous section and summarized 
by Eq.|2ni ARPES directly probes the one-particle spec- 
tral function A(k,uj). However, in extracting quantita- 
tive information from the experiment, not only the effect 
of the matrix element term /o(k, t/. A) has to be taken 
into account, but also the finite experimental resolution 
and the extrinsic continuous background due to the sec- 
ondaries (those electrons which escape from the solid af- 
ter having suffered inelastic scattering events and, there- 
fore, with a reduced Ekin)- The latter two effects may be 
explicitly accounted for by considering a more realistic 
expression for the photocurrent I(k,uj): 



dudkIo{k,i^,A)f{i:;j)A{k,tIj)R{uj~Cj)Q(k-k)+B (27) 



which consists of the convolution of Eq.jSOl with energy 
(R) and momentum (Q) resolution functions [R is typi- 
cally a Gaussian, Q may be more complicated] , and of the 
background correction B. Of the several possible forms 
for the background function B (171 . two are more fre- 
quently used: (i) the step-edge background (with three 
parameters for height, energy position, and width of the 
step-edge), which reproduces the background observed 
all the way to Ep in an unoccupied region of momentum 
space; (ii) the Shirley background Bsh{^)oc duj'P{uj'), 
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which allows to extract from the measured photocurrent 
/(w) = P{uj) + cshBsh{i^) the contribution P(w) of the 
unscattered electrons (with only the parameter csh [z3)- 
Let us now very briefly illustrate the effect of the 
matrix element term /o(k, i/, A) c>c jAf^jp, which is re- 
sponsible for the dependence of the photoemission data 
on photon energy and experimental geometry, and may 
even result in complete suppression of the intensity 
[sol IsiL Is^ . l83j |. By using the commutation relation 
hp/m ■ 



k |2 



-i[x, iJ], we can write |My - 



k\|2 



where e is a unit vector along the polarization direction 
of the vector potential A. As in Fig. [7^, let us consider 
photoemission from a dx2_y2 orbital, with the detector lo- 
cated in the mirror plane (when the detector is out of the 
mirror plane, the problem is more complicated because of 
the lack of an overall well defined even/odd symmetry). 
In order to have non vanishing photoemission intensity, 
the whole integrand in the overlap integral must be an 
even function under reflection with respect to the mirror 
plane. Because odd parity final states would be zero ev- 
erywhere on the mirror plane and therefore also at the 
detector, the final state wavefunction itself must be 
even. In particular, at the detector the photoelectron is 
described by an even parity plane-wave state e**"" with 
momentum in the mirror plane and fronts orthogonal to 
it [S^l- In turn, this implies that {e-K)\<j)f) must be even. 
In the case depicted in Fig.[7^ where is also even, the 
photoemission process is symmetry allowed for A even or 
in-plane (i.e., £p-x depends only on in-plane coordinates 
and is therefore even under reflection with respect to the 
plane) and forbidden for A odd or normal to the mirror 
plane (i.e., £s'X is odd as it depends on normal-to-the- 
plane coordinates). For a generic initial state of either 
even or odd symmetry with respect to the mirror plane, 
the polarization conditions resulting in an overall even 
matrix element can be summarized as: 



9" even 
hY odd 



A odd 



(28) 



In order to discuss the photon energy dependence, from 
Eq.^land by considering a plane wave e'*"" for the pho- 
toelectron at the detector, one may more conveniently 



write \M 



k 12 



(X 



•k)(0! 



,k|gikr\ 



The overlap integral, 



as sketched in Fig.[7|3, strongly depends on the details 
of the initial state wavefunction (peak position of the 
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FIG. 8: (Color). Beamline equipped with a plane grating 
monochromator and a 2D position-sensitive electron analyzer. 



> 
E 

3 



^ l(k,iu=.^ 



0.35 » [1« , 0.45 0.50 0.55 



K-0 «6A 
►,.1 14 evA 



-300 H 
0-25 



FIG. 9; (Color). Energy (lu) versus momentum (ky) image 
plot of the photoemission intensity from Bi2Sr2CaCu208+i 
along (0,0)-(7r,7r). This fc-space cut was taken across the Fermi 
surface (see sketch of the 2D Brillouin zone) and allows a di- 
rect visualization of the photohole spectral function A(k,ij) 
(weighted by Fermi distribution and matrix elements): the 
quasiparticle dispersion can be clearly followed up to i?F, as 
emphasized by the white circles. Energy scans at constant 
momentum (right) and momentum scans at constant energy 
(top) define energy distribution curves (EDCs) and momen- 
tum distribution curves (MDCs), respectively (from Ref. l86fl . 



radial part and oscillating character of it), and on the 
wavelength of the outgoing plane wave. Upon increas- 
ing the photon energy, both Ekin and k increase, and 
Af^j changes in a non-necessarily monotonic fashion (see 
Fig-Et; for the Cu 3d and the O 2p atomic case). In fact, 
the photoionization cross section is usually characterized 
by one minimum in free atoms, the so-called Cooper min- 
imum !84l , and a series of them in solids 183 ■ 



VI. STATE-OF-THE-ART PHOTOEMISSION 

The configuration of a generic angle-resolved photoe- 
mission beamline is shown in Fig.|S| A beam of white 
radiation is produced in a wiggler or an undulator (these 
so-called 'insertion devices' are the straight sections of 
the electron storage ring where radiation is produced); 
the light is then monochromatized at the desired photon 
energy by a grating monochromator, and is focused on 
the sample. Alternatively, a gas-discharge lamp can be 
used as a radiation source (once properly monochroma- 
tized, to avoid complications due to the presence of differ- 
ent satellites and refocused to a small spot size, essential 
for high angular resolution). However, synchrotron radi- 
ation offers important advantages: it covers a wide spec- 
tral range (from the visible to the X-ray region) with an 
intense and highly polarized continuous spectrum, while 
a discharge lamp provides only a few resonance lines at 
discrete energies. Photoemitted electrons are then col- 
lected by the analyzer, where kinetic energy and emission 
angle are determined (the whole system is in ultra-high 
vacuum at pressures lower than SxlO"""^^ torr). 

A conventional hemispherical analyzer consists of a 
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FIG. 10: Kinetic energy dependence of the 'universal' mean 
free path for excited electrons in solids (from Ref. ls^ . 



multi-element electrostatic input lens, a hemispherical 
deflector with entrance and exit slits, and an electron 
detector (i.e., a channeltron or a multi-channel detec- 
tor). The heart of the analyzer is the deflector which 
consists of two concentric hemispheres (of radius Ri 
and i?2)- These are kept at a potential difference AV, 
so that only those electrons reaching the entrance slit 
with kinetic energy within a narrow range centered at 
Epass = eAV/{Ri/R2 — R2/R1) wiU pass through this 
hemispherical capacitor, thus reaching the exit slit and 
then the detector. This way it is possible to measure 
the kinetic energy of the photoelectrons with an energy 
resolution given by AEa = Epass {w /Rq + c?l^ , where 
Ra — (^i+i?2)/2, w is the width of the entrance slit, 
and a is the acceptance angle. The role of the electro- 
static lens is that of decelerating and focusing the pho- 
toelectrons onto the entrance slit. By scanning the lens 
retarding potential one can effectively record the photoe- 
mission intensity versus the photoelectron kinetic energy. 
One of the innovative characteristics of the state-of-the- 
art analyzer is the two-dimensional position-sensitive de- 
tector consisting of two micro-channel plates and a phos- 
phor plate in series, followed by a CCD camera. In this 
case, no exit slit is required: the electrons, spread apart 
along the Y axis of the detector (Fig.lSJ as a function of 
their kinetic energy due to the travel through the hemi- 
spherical capacitor, are detected simultaneously [in other 
words, a range of electron energies is dispersed over one 
dimension of the detector and can be measured in paral- 
lel; scanning the lens voltage is in principle no longer 
necessary, at least for narrow energy windows (a few 
percent of i?j,ass)]- Furthermore, contrary to a conven- 
tional electron spectrometer in which the momentum in- 
formation is averaged over all the photoelectrons within 
the acceptance angle (typically ±1°), state-of-the-art 2D 
position-sensitive electron analyzers can be operated in 
angle-resolved mode, which provides energy-momentum 
information not only at a single fc-point but along an ex- 
tended cut in fc-space. In particular, the photoelectrons 
within an angular window of ^ 14° along the direction 



deflned by the analyzer entrance slit are focused on dif- 
ferent X positions on the detector (Fig.|SJ). It is thus 
possible to measure multiple energy distribution curves 
simultaneously for different photoelectron angles, obtain- 
ing a 2D snapshot of energy versus momentum (Fig.lSJ. 

State-of-the-art spectrometers typically allow for en- 
ergy and angular resolutions of approximately a few meV 
and 0.2°, respectively. Taking as example the transition 
metal oxides and in particular the cuprate superconduc- 
tors (for which 27r/a~ 1.6 A"-*^), one can see from Eq.[7| 
that 0.2° corresponds to ^0.5% of the Brillouin zone size, 
for the 21.2 eV photons of the Hela line typically used 
in ARPES systems equipped with a gas-discharge lamp. 
In the case of a beamline, to estimate the total energy 
resolution one has to take into account also AEm of the 
monochromator, which can be adjusted with entrance 
and exit slits (the ultimate resolution a monochromator 
can deliver is given by its resolving power R = E / AE^ ; it 
can be as good as 1-2 meV for 20 eV photons but wors- 
ens upon increasing the photon energy). To maximize 
the signal intensity at the desired total AE, monochro- 
mator and analyzer should be operated at comparable 
resolutions. As for the momentum resolution Aky, note 
that from Eq.[7| and neglecting the contribution due to 
the flnite energy resolution one can write: 

Ak|| ~ y2^^J^^~7?l2 . coszJ-Ai? (29) 

where Ai3 is the finite acceptance angle of the electron 
analyzer. From Eq.[211it is clear that the momentum res- 
olution is better at lower photon energy (i.e., lower Ekm), 
and larger polar angles 1!} (one can effectively improve the 
momentum resolution by extending the measurements to 
momenta outside the first Brillouin zone). 

Because at lower photon energies it is possible to 
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FIG. 11: High energy angle-integrated resonance photoemis- 
sion data from Ce compounds at r = 20 K (from Ref . l88fl . 
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FIG. 12: (Color). ARPES spectra and corresponding intensity plot from Sr2Ru04 along (a) F-M, and (b) M-X. (c) Measured 
and (d) calculated [89| Fermi surface. All data were taken at 10 K on a Sr2Ru04 single crystal cleaved at 180 K (from Ref. looh . 



achieve higher energy and momentum resolution, most of 
the ARPES experiments are performed in the ultraviolet 
(typically for hv < 100 eV). An additional advantage is 
that at low photon energies one can disregard the photon 
momentum k/j^ = 27r/A in Eq.[71 as for 100 eV photons 
the momentum is 0.05 A^^ (only 3% of the Brillouin zone 
size, by taking again the cuprates as an example) , and at 
21.2 eV (Hela) it is only 0.008 A^^ (0.5% of the zone). 
If on the contrary the photon momentum is not negligi- 
ble, the photoemission process does not involve vertical 
transitions and k must be explicitly taken into account in 
Eq.[71 For example, for 1487 eV photons (the Al hue 
commonly used in X-ray photoemission) kft^~0.76A~^, 
which corresponds to 50% of the zone size. 

A major drawback of working at low photon energies is 
the extreme surface sensitivity. As shown in Fig.^1 the 
mean free path for unscattered photoelectrons is charac- 
terized by a minimum of approximately 5 A at 20-100 eV 
kinetic energies ^8^] , which are typical values in ARPES 
experiments. This means that a considerable fraction 
of the total photoemission intensity will be representa- 
tive of the topmost surface layer, especially on systems 
characterized by a large structural/electronic anisotropy. 
Therefore, ARPES experiments have to be performed 
on atomically clean and well-ordered systems, which im- 
plies that atomically fresh and flat surfaces have to be 
'prepared' immediately prior to the experiment in ultra- 
high vacuum conditions (typically at pressures lower than 
5x10"^^ torr). Even then, however, because of the lower 
atomic coordination at the surface, the coexistence of 
bulk and surface electronic states, and the possible occur- 
rence of chemical and/or structural surface instabilities, 
photoemission data may not always be representative of 
the intrinsic bulk electronic structure. In order address 
with this issue, great care has to be taken also over the 
structural and chemical characterization of the sample 
surface, which can be done independently by low-energy 
electron diffraction (LEED) and core-level X-ray photoe- 
mission spectroscopy (XPS), respectively (either prior to 
or during the ARPES experiments). In this regard it has 
to be emphasized that, although the ultimate resolutions 



are not as good as in the UV regime, the sensitivity to 
bulk over surface electronic states can be enhanced (see 
Fig. llO|) by performing the ARPES experiments in the 
soft X-ray regime (500-1500 eV). The significance of this 
approach is well exemplified by recent angle-integrated 
resonance photoemission experiments performed on Ce 
compounds 88] . These Kondo systems are characterized 
by a very different degree of hybridization between the 4/ 
electronic states and other valence bands: the hybridiza- 
tion is stronger the larger the Kondo temperature Tr-. 
However, although CeRu2Si2 and CeRu2 are character- 
ized by very different Tr- (approximately 22 and 1000 
K, respectively), earlier photoemission studies reported 
similar spectra for the Ce 4/ electronic states. By per- 
forming angle-integrated high resolution photoemission 
experiments at the 3d-4/ {hv~880 eV, AS ~ 100 meV) 
and Ad-Af (ft,zy~120 eV, Ai?~50 meV) resonances (see 
Fig. lll|) . it was observed that, while the spectra for the 
two compounds are indeed qualitatively similar at 120 eV 
photon energy, they are remarkably different at 880 eV. 
As the photoelectron mean free path increases from ap- 
proximately 5 to almost 20 A upon increasing the photon 
energy from 120 to 880 eV (Fig. lTUIl . it was concluded that 
the 4c?-4/ spectra mainly reflect the surface 4/ electronic 
states. These are different from those of the bulk and are 
not representative of the intrinsic electronic properties of 
the two compounds, which are more directly probed at 
880 eV: the 3c?-4/ spectra show a prominent structure 
corresponding to the tail of a Kondo peak in CeRu2Si2, 
and a broader feature reflecting the more itinerant char- 
acter of the 4/ electrons in CeRu2 [s^. 

In the following, we will move on to the review of recent 
ARPES results from several materials, such as Sr2Ru04, 
2i/-NbSe2, Be(OOOl), and Mo(llO). These examples wifl 
be used to illustrate the capability of this technique and 
some of the speciflc issues that one can investigate in de- 
tail by ARPES. In particular, these test cases will demon- 
strate that, by taking full advantage of the momentum 
and energy resolution as well as of the photon energy 
range nowadays available, state-of-the-art ARPES is a 
unique tool for momentum space microscopy. 
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A. Sr2Ru04: Bands and Fermi Surface 

To illustrate how one can study electronic bands and 
Fermi surfaces by ARPES, and how critical the improve- 
ment in resolution has been in this regard, the novel 
superconductor Sr2Ru04 is a particularly good exam- 
ple. Its low-energy electronic structure, as predicted 
by band-structure calculations is characterized by three 
bands crossing the chemical potential [oj- These 
define a complex Fermi surface comprised of two elec- 
tron pockets and one hole pocket fFig.ll2tl'l. which have 
been clearly observed in de Haas- van Alphen experiments 
Iq^ I . On the other hand, early photoemission mea- 
surements suggested a different topology 95, 96, 97], 
which generated a certain degree of controversy in the 
field j98| . This issue was conclusively resolved only by 
taking advantage of the high energy and momentum res- 
olution of the 'new generation' of ARPES data: it was 
then recognized that a surface reconstruction and, 
in turn, the detection of several direct and folded sur- 
face bands were responsible for the conflicting interpre- 
tations 113, Ens Em Eo3- Fig.[Tli,b show high res- 
olution ARPES data {AE = 14 meV, Afc = 1.5% of 
the zone edge) taken at 10 K with 28 eV photons on a 
Sr2Ru04 single crystal cleaved at 180 K (for Sr2Ru04, 
as recently discovered, high-temperature cleaving sup- 
presses the reconstructed-surface contributions to the 
photoemission signal and allows one to isolate the bulk 
electronic structure Many well defined quasiparti- 

cle peaks disperse towards the Fermi energy and disap- 
pear upon crossing Ep- A Fermi energy intensity map 
(Fig.Et) can then be obtained by integrating the spec- 
tra over a narrow energy window about Ep (±10 meV). 
As the spectral function (multiplied by the Fermi func- 
tion) reaches its maximum at Ep when a band crosses the 
Fermi energy, the Fermi surface is identified by the local 
maxima of the intensity map. Following this method, the 
three sheets of Fermi surface are clearly resolved and are 
in excellent agreement with the theoretical calculations 
(Fig.CH). 



B. 2_ff-NbSe2: Superconducting Gap 




FIG. 13: (Color), (a) 2//-NbSe2 ARPES spectra (measured 
at 20 K with 21.2 eV photons), (b) corresponding image plot, 
and (c) band structure calculations along T-K (from Ref. 105) . 



of theoretical calculations (as shown in Fig.^J where 
ARPES spectra and band structure calculations are com- 
pared for the T-K high symmetry direction). As for the 
low temperature charge-density wave phase, despite the 
intense effort no agrement has been reached yet o n the 
driving force responsible for the transition [33, Il04| . 

Owing to the great improvement in energy and mo- 
mentum resolution, it has now become possible to study 
by ARPES also the momentum and temperature depen- 
dence of the superconducting gap on low-Tj, materials 
(until recently, experiments of this kind could been per- 
formed only for the much larger d-wave gap of the high- 
Tc superconductors 0|). The data presented in Fig. [HI 
which are one of the most impressive examples of com- 
bined high energy and momentum resolution in ARPES 
experiments on solid samples (i.e., AE = 2.5 meV and 
Afc = 0.2°), provide direct evidence for Fermi surf ace 
sheet-dependent superconductivity in 2iJ-NbSe2 '106*1 . A 
superconducting gap of about 1 meV was successfully de- 
tected along two of the normal-state Fermi surface sheets, 
but not along the third one. In fact, the opening of the 
gap is directly evidenced in Fig. [Tib andlTlh by the shift 



2iJ-NbSe2 is an interesting quasi two-dimensional sys- 
tem exhibiting a charge-density wave phase transition 
at approximately 33 K, and a phonon-mediated super- 
conducting phase transitio n at 7.2 K. As indicated by 
band structure calculations 'l03], the valence-band elec- 
tronic structure is characterized by a manifold of disper- 
sive bands in a 6 eV range below the Fermi energy. At 
low energy, three dispersive bands are expected to cross 
the chemical potential and define three sheets of Fermi 
surface in the hexagonal Brillouin zone. Both the band 
manyfold and the Fermi surface topology have been stud- 
ied in great detail by ARPES; exception made for a weak 
energy renormalization, the normal-state experimental 
data are in extremely good agreement with the results 
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FIG. 14: (Color). Normal (10 K, red) and superconducting 
state (5.3 K, blue) ARPES spectra from 2//-NbSe2, measured 
at fc-points belonging to the three different sheets of Fermi 
surface (see insets). The value of the superconducting gap 
used to fit the data is indicated in each panel (from Ref. ilO&l . 
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to high binding energies of the 5.3 K spectra leading-edge 
midpoint (which is instead located at Ep at 10 K, as ex- 
pected for a metal), and by the simultaneous appearance 
of a peak below Ep (which reflects the piling up of the 
density of states due to the gap opening). This behavior 
is absent for the inner Fermi surface pocket fFig.ll4b^. 

C. Self Energy and Collective Modes 

As discussed in Sec. II VI the introduction of the elec- 
tron self energy E(k,a;) = E'(k,a;)-fiI]"(k, w) is a power- 
ful way to account for many-body correlations in solids. 
Its real and imaginary parts correspond, respectively, to 
the energy renormalization with respect to the bare band 
energy and to the finite lifetime of the quasiparticles in 
the interacting system. Owing to the energy and momen- 
tum resolution nowadays achievable, both components of 
the self energy can be in principle estimated very accu- 
rately from the analysis of the ARPES intensity in terms 
of energy distribution curves (EDCs) and/or momentum 
distribution curves (MDCs) , which is one of the aspects 
that make ARPES such a powerful tool for the investiga- 
tion of complex materials. In some cases the MDC anal- 
ysis may be more effective than the analysis of the EDCs 
in extracting information on the self energy. In fact, 
EDCs are typically characterized by a complex lineshape 
(Fig.O because of the nontrivial u! dependence of the 
self energy, the presence of additional background, and 
the low-energy cutoff due to the Fermi function. Further- 
more, as evidenced by the generic expression for the spec- 
tral function A(k, oj) in Eg . 1221 the EDC peak position 
is determined by I]'(k, w) as well as S]"(k, w), because 
both terms are strongly energy dependent. On the other 
hand, if the self energy is independent of k normal to 
the Fermi surface (and the matrix elements are a slowly- 
varying function of fc), then the corresponding MDCs 
are Lorentzians centered at k = kp + [oJ — (to)]/ Vp with 
FWHM given by 2'E"{uj)/vp, where Vp is the bare Fermi 
velocity normal to the Fermi surface [this is obtained by 
approximating — Vpik — kp) in Ea.l22|. Lorentzian 
lineshapes were indeed observed for the MDCs (Fig.|5Jl. 

As an example of this kind of analysis we will briefly 
discuss the case of electron-phonon coupling on metallic 
surfaces, for which the established t heoretical formalism 
can be applied very effectively [M [lOl [lOl [lOl [lT(l 
llllUll2l |. The electron-phonon interaction involving sur- 
face phonons and the F-surface state on the Be(OOOl) 
surface was investigated by two g roups, and qualita - 
tively similar conclusion were 

drawn [Tol Hoi Efl 
Fig. 115b shows results for the Be(OOOl) surface state along 
the TM direction of the surface Brillouin z one; a feature 
is seen dispersing towards the Fermi level |l09l | . Close 
to Ep the spectral function exhibits a complex structure 
characterized by a broad hump and a sharp peak, with 
the latter being confined to within an energy range given 
by the typical bandwidth LOph of the surface phonons. 
This behavior corresponds to a 'two-branch' splitting of 



the neav-Ep dispersion, with a transfer of spectral weight 
between the two branches as a function of binding en- 
ergy. While the high-energy dispersion is representative 
of the bare quasiparticles, at low energy the dispersion 
is renormalizcd by the electron-phonon interaction (this 
behavior is shown, for a similar electron-phonon coupled 
system, in the inset of Fig.ll5b'). In other words, the 
weaker dispersion observed at energies smaller than ujph 
describes dressed quasiparticles with an effective mass 
enhanced by a factor of (1-l- A), w here A is the electron- 
phonon coupling parameter |114| . The latter can also 
be estimated from the ratio of renormalizcd (vj^) and 
bare (v^) quasiparticle velocities, according to the rela- 
tion Vj^ = /i^^9£k/9k= (1+A)^^v|^. This way, for the data 
presented in Fig.ll5b. the value A = 1.18 was obtained (al- 
ternatively A can also be estimated fro m th e temperature 
dependence of the linewidth near Ep |l07j ). 

A similar example of electron-phonon coupled system 
is the surface of Mo(llO) |113| . In this case, the real and 
imaginary part of the self energy shown in Fig. 115b were 
obtained directly from the EDC analysis: E" corresponds 
to the EDC width and E' to the difference between the 
observed quasiparticle dispersion and a straight line ap- 
proximating the dispersion of the non-interacting system 
fFig.ll5b. inset). The step-like change at 30 meV in E" 
is interpreted as the phonon contribution (dashed line) 
and the parabolic part at higher energies is attributed to 
electron-electron interactions. The phonon contribution 
to the real part of the self energy is calculated from the 
Kramers-Kronig relations (see Sec. lIV|l and agrees well 
with the data (dotted line). As an additional confirma- 
tion of the electron-phonon description, it was noted that 
the temperature dependence of the s cattering rate is well 
reproduced by the calculations |113| . 




ly^'i: mo hoq iw 2i'\). Cr 4(j"j -hj iui lot « 



l.nerxjv '.mpM'/ E-E|- ImcV) 

FIG. 15: (a) ARPES spectra for the Be(OOOl) surface state 
(from R,ef. lT09|) . (b) Self energy estimated from the Mo(llO) 
surface state ARPES spectra, and corresponding quasiparticle 
dispersion (inset). Calculated electron-phonon contributions 
to the real and imaginary part of E(k, lj) are indicated by 
dotted and dashed lines, respectively (the latter was off set b y 
26meV to account for impurity scattering). From Ref. lllA 



14 



[1] H. Hertz, Ann. Phys. 17, 983 (1887). 
[2] A. Einstein, Ann. Physik 31, 132 (1905). 
[3] R. Z. Bachrach, Synchrotron Radiation Research, Ad- 
vances in Surface and Interface Science, vol. 1 (Plenum 
Press, New York, 1992). 
[4] J. Braun, Rep. Prog. Phys. 59, 1267 (1996). 
[5] C. R. Brundle and A. D. Baker, Electron Spectroscopy: 
Theory, Techniques, and Applications, vol. 1 (Academic 
Press, New York, 1977). 
[6] C. R. Brundle and A. D. Baker, Electron Spectroscopy: 
Theory, Techniques, and Applications, vol. 2 (Academic 
Press, New York, 1978). 
[7] M. Cardona and L. Ley, Photoemission in Solids, vol. 1 

( Springer- Verlag, Berlin, 1978). 
[8] T. A. Carlson, Photoelectron and Auger Spectroscopy 

(Plenum Press, New York, 1975). 
[9] R. Courths and S. Hiifner, Phys. Rep. 112, 53 (1984). 
[10] A. Damascelli, D. H. Lu, and Z.-X. Shen, J. Electron 

Spectr. Relat. Phenorn. 117-118, 165 (2001). 
[11] A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. 

Phys. 75, 473 (2003). 
[12] D. E. Eastman, in Techniques of Metal Research, edited 
by E. Passaglia (Interscience Publisher, New York, 
1972), vol. VI, part 1. 
[13] B. Feuerbacher and R. F. Willis, J. Phys. Solid State 

Phys. 9, 169 (1976). 
[14] B. Feuerbacher, B. Fitton, and R. F. Willis, Photoe- 
mission on the Electronic Properties of Surfaces (John 
Wiley & Sons, New York, 1978). 
[15] M. Grioni, J. Electron Spectr. Relat. Phenom. 117-118 

(2001), special issue on Strongly Correlated Systems. 
[16] F. J. Himpsel, Adv. Phys. 32, 1 (1983). 
[17] S. Hiifner, Photoelectron Spectroscopy (Springer- Verlag, 
Berlin, 1995). 

[18] J. E. Inglesfield and B. W. Holland, in The Chemical 

Physics of Solid Surfaces and Heterogeneous Catalysis, 
edited by D. A. King and D. P. Woodruff (Elsevier Sci- 
entific Publishing Co., Amsterdam, 1981), vol. 1. 
[19] S. D. Kevan, Angle Resolved Photoemission- Theory and 
Current Applications (Elsevier Science, Amsterdam, 
1992). 

[20] R. C. G. Leckey, Appl. Surf. Sci. 13, 125 (1982). 

[21] L. Ley and M. Cardona, Photoemission in Solids, vol. II 
(Springer- Verlag, Berlin, 1979). 

[22] I. Lindau and W. E. Spicer, in Synchrotron Radiation 
Research, edited by H. Winick and S. Doniach (Plenum 
Press, New York, 1980). 

[23] D. W. Lynch and C. G. Olson, Photoemission Studies 
of High- Temperature Superconductors (Cambridge Uni- 
versity Press, Cambridge, 1999). 

[24] G. D. Mahan, in Electron and Ion Spectroscopy 
of Solids, edited by L. Fiermans, J. Vennik, and 
W. Dekeyser (Plenum Press, New York, 1978). 

[25] G. Margaritondo and J. H. Weaver, in Methods of Ex- 
perimental Physics: Surfaces, edited by M. G. Legally 
and R. L. Park (Academic Press, New York, 1983). 

[26] V. V. Nemoshkalenko and V. G. Aleshin, Electron Spec- 
troscopy of Crystals (Plenum Press, New York, 1979). 

[27] E. W. Plummer and W. Eberhardt, in Advances in 
Chemical Physics, edited by I. Prigogine and S. A. Rice 
(John Wiley & Sons, New York, 1982), vol. XLIV. 



[28] Z.-X. Shen and D. S. Dessau, Phys. Rep. 253, 1 (1995). 

[29] N. V. Smith, Crit. Rev. Solid State Sci. 2, 45 (1971). 

[30] N. V. Smith and F. J. Himpsel, in Handbook on 
Synchrotron Radiation, edited by E.-E. Koch (North- 
Holland, Amsterdam, 1983), vol. lb. 

[31] K. E. Smith and S. D. Kevan, Prog. Sohd State Chem. 
21, 49 (1991). 

[32] G. Wendin, Breakdown of the One-Electron Pictures 
in Photoelectron Spectroscopy (Springer- Verlag, Berlin, 
1981). 

[33] G. Wertheim, in Electron and Ion Spectroscopy of Solids, 

edited by L. Fiermans, J. Vennik, and W. Dekeyser 

(Plermm Press, New York, 1978). 
[34] R. H. Williams, G. P. Srivastava, and I. T. McGovern, 

Rep. Prog. Phys. 43, 1357 (1980). 
[35] M. B. J. Meinders, Ph.D. Thesis (1994), University of 

Groningen, The Netherlands. 
[36] G. D. Mahan, Phys. Rev. B 2, 4334 (1970). 
[37] T. Pillo, Ph.D. Thesis (1999), University of Freiburg, 

Switzerland. 

[38] V. N. Strocov, H. I. Starnberg, P. O. Nilsson, H. E. 
Brauer, and L. J. HoUeboom, Phys. Rev. Lett. 79, 467 
(1997). 

[39] V. N. Strocov, R. Claessen, G. Nicolay, S. Hiifner, 
A. Kimura, A. Harasawa, S. Shin, A. Kakizaki, P. O. 
Nilsson, H. I. Starnberg, et al., Phys. Rev. Lett. 81, 
4943 (1998). 

[40] E. D. Hansen, T. Miller, and T.-C. Chiang, Phys. Rev. 

Lett. 80, 1766 (1998). 
[41] N. V. Smith, P. Thiry, and Y. Petroff, Phys. Rev. B 47, 

15476 (1993). 

[42] T. Miller, W. E. McMahon, and T.-C. Chiang, Phys. 

Rev. Lett. 77, 1167 (1996). 
[43] E. D. Hansen, T. Miller, and T.-C. Chiang, Phys. Rev. 

B 55, 1871 (1997). 
[44] E. D. Hansen, T. Miller, and T.-C. Chiang, Phys. Rev. 

Lett. 78, 2807 (1997). 
[45] K. Mitchell, Proc. Roy. Soc. London A 146, 442 (1934). 
[46] R. E. B. Makinson, Phys. Rev. 75, 1908 (1949). 
[47] M. J. Buckingham, Phys. Rev. 80, 704 (1950). 
[48] W. L. Schaich and N. W. Ashcroft, Phys. Rev. B 3, 

2452 (1971). 

[49] P. J. Feibelman and D. E. Eastman, Phys. Rev. B 10, 

4932 (1974). 
[50] J. B. Pendry, Surf. Sci. 57, 679 (1976). 
[51] J. B. Pendry, J. Phys. Solid State Phys. 8, 2413 (1975). 
[52] A. Liebsch, Phys. Rev. B 13, 544 (1976). 

[53] A. Liebsch, in Electron and Ion Spectroscopy of Solids, 
edited by L. Fiermans, J. Vennik, and W. Dekeyser 
(Plenum Press, New York, 1978). 

[54] M. Lindroos and A. Bansil, Phys. Rev. Lett. 75, 1182 

(1995) . 

[55] M. Lindroos and A. Bansil, Phys. Rev. Lett. 77, 2985 

(1996) . 

[56] A. Bansil and M. Lindroos, J. Phys. Chem. Solids 56, 

1855 (1995). 

[57] A. Bansil and M. Lindroos, J. Phys. Chem. Solids 59, 
1879 (1998). 

[58] A. Bansil and M. Lindroos, Phys. Rev. Lett. 83, 5154 

(1999). 

[59] H. Y. Fan, Phys. Rev. 68, 43 (1945). 



15 



[60] C. N. Berglund and W. E. Spicer, Phys. Rev. 136, 

A1030 (1964). [90 

[61] G. A. Sawatzky, Nature 342, 480 (1989). 

[62] J. W. Gadzuk and M. Sunjic, Phys. Rev. B 12, 524 

(1975). [91 

[63] K. Siegbahn, K. Siegbahn, C. Nordling, G. Johans- [92; 
son, J. Hedman, P. F. Heden, K. Hamrin, U. I. Gelius, [93 
T. Bergmark, L. O. Wcrrnc, ot al., ESCA Applied to 
Free Molecules (North Holland, Amsterdam, 1969). 

[64] A. A. Abrikosov, L. P. Gor'kov, and I. E. Dzyaloshin- [94 
skii. Quantum Field Theoretical Methods in Statistical 
Physics (Pergamon Press, Oxford, 1965). 

[65] L. Hedin and S. Lundqvist, in Solid State Physics: [95 
Advances in Research and Applications, edited by 
H. Ehrcnreich, F. Seitz, and D. Turnbull (Academic, 
New York, 1969), vol. 23. [96; 

[66] A. L. Fetter and J. D. Walccka, Quantum Theory 
of Many-Particle Systems (McGraw-Hill, New York, 
1971). [97 

[67] G. D. Mahan, Many-Particle Physics (Plenum Press, 
New York, 1981). 

[68] E. N. Economou, Green's Functions in Quantum [98 
Physics, vol. 7 of Springer Series in Solid State Science 
(Springer- Verlag, Berlin, 1983). [99; 

[69] G. Rickayzcn, Green's Functions and Gondensed Matter 
in Techniques of Physics, vol. Vol. 7 (Academic Press, 
London, 1991). [lOO; 

[70] M. R. Norman, H. Ding, H. FrotwcU, M. Randeria, and 
J. C. Campuzano, Phys. Rev. B 60, 7585 (1999). 

[71] L. D. Landau, Sov. Phys. JETP 3, 920 (1956). 

[72] L. D. Landau, Sov. Phys. JETP 5, 101 (1957). [101 

[73] L. D. Landau, Sov. Phys. JETP 8, 70 (1959). 

[74] D. Pines and P. Noziores, The Theory of Quantum Liq- [102 
uids, vol. 1 (Benjamin, New York, 1966). 

[75] J. M. Luttinger, Phys. Rev. 121, 942 (1961). 

[76] C. Hodges, H. Smith, and J. W. Wilkins, Phys. Rev. B [103 

4, 302 (1971). 

[77] J. J. Yeh and I. Lindau, Atom. Data Nucl. Data Tabl. 

32, 1 (1985). [104 

[78] P. Nozieres, Theory of Interacting Fermi Systems (Ben- 
jamin, New York, 1964). 

[79] D. A. Shirley, Phys. Rev. B 5, 4709 (1972). [105 

[80] G. W. Gobeh, F. G. Allen, and E. O. Kane, Phys. Rev. [106^ 
Lett. 12, 94 (1964). 

[81] E. Dictz, H. Becker, and U. Gerhardt, Phys. Rev. Lett. [107 
36, 1397 (1976). 

[82] J. Hermanson, Solid State Commun. 22, 9 (1977). [108 
[83] W. Eberhardt and F. J. Himpsel, Phys. Rev. B 21, 5572 

(1980). [109 
[84] J. W. Cooper, Phys. Rev. 128, 681 (1962). 

[85] S. L. Molodtsov, S. V. Halilov, V. D. P. Serve- [lio; 
dio, W. Schneider, S. Danzenbacher, J. J. Hinarejos, 
M. Richter, and C. Laubschat, Phys. Rev. Lett. 85, 4184 [111 
(2000). 

[86] T. Valla, A. V. Fedorov, P. D. Johnson, B. O. Wells, [112; 

5. L. Hulbert, Q. Li, G. D. Gu, and N. Koshizuka, Sci- 
ence 285, 2110 (1999). [113 

[87] M. P. Seah and W. A. Dench, Surf. Interface Anal. 1,2 

(1979). [114; 

[88] A. Sekiyama, T. Iwasaki, K. Matsuda, Y. Saitoh, 
Y. Onuki, and S. Suga, Nature 403, 396 (2000). 

[89] I. I. Mazin and D. J. Singh, Phys. Rev. Lett. 79, 733 



(1997). 

A. Damascelli, D. H. Lu, K. M. Shen, N. P. Armitage, 
F. Ronning, D. L. Feng, C. Kim, Z.-X. Shen, T. Kirrmra, 
Y. Tokura, et al., Phys. Rev. Lett. 85, 5194 (2000). 
D. J. Singh, Phys. Rev. B 52, 1358 (1995). 
T. Oguchi, Phys. Rev. B 51, 1385 (1995). 
A. P. Mackenzie, S. R. Julian, A. J. Diver, G. J. Mcmul- 
lan, M. P. Ray, G. G. Lonzarich, Y. Maeno, S. Nishizaki, 
and T. Fujita, Phys. Rev. Lett. 76, 3786 (1996). 

C. Bergemann, S. R. Julian, A. P. Mackenzie, 
S. NishiZaki, and Y. Maeno, Phys. Rev. Lett. 84, 2662 
(2000). 

T. Yokoya, A. Chainani, T. Takahashi, H. Ding, J. C. 
Campuzano, H. Katayama-Yoshida, M. Kasai, and 
Y. Tokura, Phys. Rev. B 54, 13311 (1996). 
T. Yokoya, A. Chainani, T. Takahashi, H. Kayatama- 
Yoshida, M. Kasai, and Y. Tokura, Phys. Rev. Lett. 76, 
3009 (1996). 

D. H. Lu, M. Schmidt, T. R. Cummins, S. Schuppler, 
F. Lichtenberg, and J. G. Bednorz, Phys. Rev. Lett. 76, 

4845 (1996). 

A. V. Puchkov, Z.-X. Shen, T. Kimura, and Y. Tokura, 
Phys. Rev. B 58, 13322 (1998). 

R. Matzdorf, Z. Fang, Ismail, J. Zhang, T. Kimura, 
Y. Tokura, K. Terakura, and E. W. Plummer, Science 
289, 746 (2000). 

A. Damascelh, K. M. Shen, D. H. Lu, N. P. Armitage, 
F. Ronning, D. L. Feng, C. Kim, Z.-X. Shen, T. Kimura, 
Y. Tokura, et al., J. Electron Spectr. Relat. Phenom. 
114-116, 641 (2001). 

A. Damascelh, K. M. Shen, D. H. Lu, and Z.-X. Shen, 
Phys. Rev. Lett. 87, 239702 (2001). 

K. M. Shen, A. Damascelli, D. H. Lu, N. P. Armitage, 
F. Ronning, D. L. Feng, C. Kim, Z.-X. Shen, D. J. Singh, 
I. I. Mazin, et al., Phys. Rev. B 64, 180502 (2001). 
R. Corcoran, P. Meeson, Y. Onuki, P.-A. Probst, 
M. Springford, K. Takita, H. Harima, G. Guo, and 

B. Gyorffy, J. Phys. Cond. Matt. 6, 4479 (1994). 

T. Straub, T. Finteis, R. Claessen, P. Steiner, P. Hiifner, 
Blalia, C. Oglesby, and E. Bucher, Phys. Rev. Lett. 82, 
4504 (1999). 

A. Damascelli (2000), unpublished. 

T. Yokoya, T. Kiss, A. Chainani, S. Shin, M. Nohara, 

and H. Takagi. Science 294, 2518 (2001). 

T. Balasubranianian, E. Jensen, X. L. Wu, and S. L. 

Hulbert, Phys. Rev. B 57, 6866 (1998). 

P. Hofmann, Y. Q. Cai, C. Griitter, and J. H. Bilgram, 

Phys. Rev. Lett. 81, 1670 (1998). 

M. Hengsberger, D. Purdie, P. Segovia, M. Garnier, and 

Y. Baer, Phys. Rev. Lett. 83, 592 (1999). 

M. Hengsberger, R. Fresard, D. Purdie, P. Segovia, and 

Y. Baer, Phys. Rev. B 60, 10796 (1999). 

S. LaShell, E. Jensen, and T. Balasubramanian, Phys. 

Rev. B 61, 2371 (2000). 

E. Rotenberg, J. Schaefer, and S. D. Keven, Phys. Rev. 
Lett. 84, 2925 (2000). 

T. VaUa, A. V. Fedorov, P. D. Johnson, and S. L. Hul- 
bert, Phys. Rev. Lett. 83, 2085 (1999). 
N. W. Ashcroft and N. D. Mermin, Solid State Physics 
(Saunders College, Philadelphia, 1976). 



